Search for: All records

Abstract The end of supernova remnant (SNR) evolution is characterized by a so-called “radiative” stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at the ends of their lives. A key prediction of SNR evolutionary models is the formation at the onset of the radiative stage of a cold, dense shell behind the forward shock. However, searches for these shells via their neutral hydrogen emission have had limited success. We instead introduce an independent observational signal of shell formation arising from the interaction between nonthermal particles accelerated by the SNR forward shock (cosmic rays) and the dense shell. Using a semi-analytic model of particle acceleration based on state-of-the-art simulations coupled with a high-resolution hydrodynamic model of SNR evolution, we predict the nonthermal emission that arises from this interaction. We demonstrate that the onset of the radiative stage leads to nonthermal signatures from radio to gamma rays, including radio and gamma-ray brightening by nearly 2 orders of magnitude. Such a signature may be detectable with current instruments, and will be resolvable with the next generation of gamma-ray telescopes (namely, the Cherenkov Telescope Array). 

ABSTRACT Arkenstone is a new scheme that allows multiphase, stellar feedback-driven winds to be included in coarse resolution cosmological simulations. The evolution of galactic winds and their subsequent impact on the circumgalactic medium are altered by exchanges of mass, energy, momentum, and metals between their component phases. These exchanges are governed by complex, small-scale physical processes that cannot be resolved in cosmological simulations. In this second presentation paper, we describe Arkenstone’s novel cloud particle approach for modelling unresolvable cool clouds entrained in hot, fast winds. This general framework allows models of the cloud–wind interaction, derived from state-of-the-art high-resolution simulations, to be applied in a large-scale context. In this work, we adopt a cloud evolution model that captures simultaneous cloud mass loss to and gain from the ambient hot phase via turbulent mixing and radiative cooling, respectively. We demonstrate the scheme using non-cosmological idealized simulations of a galaxy with a realistic circumgalactic medium component, using the arepo code. We show that the ability of a high-specific energy wind component to perform preventative feedback may be limited by heavy loading of cool clouds coupled into it. We demonstrate that the diverging evolution of clouds of initially differing masses leads to a complex velocity field for the cool phase and a cloud mass function that varies both spatially and temporally in a non-trivial manner. These latter two phenomena can manifest in the simulation because of our choice of a Lagrangian discretization of the cloud population, in contrast to other proposed schemes. 

Abstract We present a suite of six high-resolution chemodynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback among the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation Ω =v_circ/Rvia the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with Ω, so that the bulge-dominated galaxies have higher midplane densities, lower velocity dispersions, and higher molecular gas fractions than their main-sequence counterparts. The star formation rate in five out of six galaxies is independent of Ω and is consistent with regulation by the midplane gas pressure alone. In the sixth galaxy, which has the most centrally concentrated bulge and thus the highest Ω, we reproduce dynamical suppression of the star formation efficiency in agreement with observations. This produces a transition away from pressure-regulated star formation. 

Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over $Z\prime \equiv Z/Z_{\odot }\sim 0.1-3$, gas surface density Σ_gas∼ 5–150M_⊙pc⁻², and stellar surface density Σ_star∼ 1–50M_⊙pc⁻². The range of emergent star formation rate surface density is Σ_SFR∼ 10⁻⁴–0.5M_⊙kpc⁻²yr⁻¹, and ISM total midplane pressure isP_tot/k_B= 10³–10⁶cm⁻³K, withP_totequal to the ISM weight $\mathcal{W}$. For given Σ_gasand Σ_star, we find ${\mathrm{\Sigma }}_{\mathrm{SFR}}\propto Z{\prime }^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to Σ_SFR. The thermal feedback yield varies sensitively with both weight and metallicity as ${\mathrm{\Upsilon }}_{\mathrm{th}}\propto {\mathcal{W}}^{-0.46}Z{\prime }^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence ${\mathrm{\Upsilon }}_{\mathrm{turb}+\mathrm{mag}}\propto {\mathcal{W}}^{-0.22}Z{\prime }^{-0.18}$. The reduction in Σ_SFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 

Abstract We present the Local GroupL-Band Survey, a Karl G. Jansky Very Large Array (VLA) survey producing the highest-quality 21 cm and 1–2 GHz radio continuum images to date, for the six VLA-accessible, star-forming, Local Group galaxies. Leveraging the VLA’s spectral multiplexing power, we simultaneously survey the 21 cm line at high 0.4 km s⁻¹velocity resolution, the 1–2 GHz polarized continuum, and four OH lines. For the massive spiral M31, the dwarf spiral M33, and the dwarf irregular galaxies NGC 6822, IC 10, IC 1613, and the Wolf–Lundmark–Melotte Galaxy, we use all four VLA configurations and the Green Bank Telescope to reach angular resolutions of <5″ (10–20 pc) for the 21 cm line with <10²⁰cm⁻²column density sensitivity, and even sharper views (<2″; 5–10 pc) of the continuum. Targeting these nearby galaxies (D ≲ 1 Mpc) reveals a sharp, resolved view of the atomic gas, including 21 cm absorption, and continuum emission from supernova remnants and Hiiregions. These data sets can be used to test theories of the abundance and formation of cold clouds, the driving and dissipation of interstellar turbulence, and the impact of feedback from massive stars and supernovae. Here, we describe the survey design and execution, scientific motivation, data processing, and quality assurance. We provide a first look at and publicly release the wide-field 21 cm Hidata products for M31, M33, and four dwarf irregular targets in the survey, which represent some of the highest-physical-resolution 21 cm observations of any external galaxies beyond the LMC and SMC. 

Abstract Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of V. Springel & L. Hernquist (hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially resolved simulations of the turbulent, multiphase, star-forming interstellar medium (ISM) have begun to produce new first-principles models, which when fully developed can replace traditional subgrid prescriptions. This approach has advantages of being physically motivated and predictive rather than empirically tuned, and allowing for varying environmental conditions rather than being tied to local-Universe conditions. As a prototype of this new approach, by combining calibrations from the TIGRESS numerical framework with the pressure-regulated feedback-modulated (PRFM) theory, simple formulae can be obtained for both the gas depletion time and an effective equation of state. Considering galaxies in TNG50, we compare the “native” simulation outputs with postprocessed predictions from PRFM. At TNG50 resolution, the total midplane pressure is nearly equal to the total ISM weight, indicating that galaxies in TNG50 are close to satisfying vertical equilibrium. The measured gas scale height is also close to theoretical equilibrium predictions. The slopes of the effective equations of states are similar, but with effective velocity dispersion normalization from SH03 slightly larger than that from current TIGRESS simulations. Because of this and the decrease in PRFM feedback yield at high pressure, the PRFM model predicts shorter gas depletion times than the SH03 model at high densities and redshift. Our results represent a first step toward implementing new, numerically calibrated subgrid algorithms in cosmological galaxy formation simulations. 

Measuring the properties of the cold neutral medium (CNM) in low-metallicity galaxies provides insight into heating and cooling mechanisms in early Universe-like environments. We report detections of two localized atomic neutral hydrogen (Hi) absorption features in NGC 6822, a low-metallicity (0.2 Z⊙) dwarf galaxy in the Local Group. These are the first unambiguous CNM detections in a low-metallicity dwarf galaxy outside the Magellanic Clouds. The Local Group L-Band Survey (LGLBS) enabled these detections due to its high spatial (15 pc for Hi emission) and spectral (0.4 km s−1) resolution. We introduce LGLBS and describe a custom pipeline to search for Hi absorption at high angular resolution and extract associated Hi emission. A detailed Gaussian decomposition and radiative transfer analysis of the NGC 6822 detections reveals five CNM components, with key properties: a mean spin temperature of 32±6 K, a mean CNM column density of 3.1×1020 cm−2, and CNM mass fractions of 0.33 and 0.12 for the two sightlines. Stacking non-detections does not reveal low-level signals below our median optical depth sensitivity of 0.05. One detection intercepts a star-forming region, with the Hi absorption profile encompassing the CO (2−1) emission, indicating coincident molecular gas and a depression in high-resolution Hi emission. We also analyze a nearby sightline with deep, narrow Hi self-absorption dips, where the background warm neutral medium is attenuated by intervening CNM. The association of CNM, CO, and Hα emissions suggests a close link between the colder, denser Hi phase and star formation in NGC 6822. 

Abstract We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new GPU-accelerated version of the Athena++ AMR code, and they span more than six orders of magnitude in radius, reaching scales similar to that of the black hole horizon. The key physical ingredients are radiative cooling and a phenomenological heating model. We find that the accretion flow takes the form of multiphase gas at radii less than about a kpc. The cold gas accretion includes two dynamically distinct stages: the typical disk stage in which the cold gas resides in a rotationally supported disk, and relatively rare chaotic stages (≲10% of the time) in which the cold gas inflows via chaotic streams. Though cold gas accretion dominates the time-averaged accretion rate at intermediate radii, accretion at the smallest radii is dominated by hot virialized gas at most times. The accretion rate scales with radius as M ̇ ∝ r 1 / 2 when hot gas dominates, and we obtain M ̇ ≃ 10 − 4 – 10 − 3 M ⊙ yr − 1 near the event horizon, similar to what is inferred from EHT observations. The orientation of the cold gas disk can differ significantly on different spatial scales. We propose a subgrid model for accretion in lower-resolution simulations in which the hot gas accretion rate is suppressed relative to the Bondi rate by ∼ ( r g / r Bondi ) 1 / 2 . Our results can also provide more realistic initial conditions for simulations of black hole accretion at the event horizon scale. 

Abstract Measuring the properties of the cold neutral medium (CNM) in low-metallicity galaxies provides insights into heating and cooling mechanisms in early Universe-like environments. We report detections of two localized atomic neutral hydrogen (Hi) absorption features in NGC 6822, a low-metallicity (0.2Z_⊙) dwarf galaxy in the Local Group. These are the first unambiguous CNM detections in a low-metallicity dwarf galaxy outside the Magellanic Clouds. The Local GroupL-band Survey (LGLBS) enabled these detections, due to its high spatial (15 pc for Hiemission) and spectral (0.4 km s⁻¹) resolution. We introduce LGLBS and describe a custom pipeline for searching for Hiabsorption at high angular resolution and extracting associated Hiemission. A detailed Gaussian decomposition and radiative transfer analysis of the NGC 6822 detections reveals five CNM components, with key properties: a mean spin temperature of 32 ± 6 K, a mean CNM column density of 3.1 × 10²⁰cm⁻², and CNM mass fractions of 0.33 and 0.12 for the two sightlines. Stacking nondetections does not reveal low-level signals below our median optical depth sensitivity of 0.05. One detection intercepts a star-forming region, with the Hiabsorption profile encompassing the CO (2−1) emission, indicating coincident molecular gas and a depression in high-resolution Hiemission. We also analyze a nearby sightline with deep, narrow Hiself-absorption dips, where the background warm neutral medium is attenuated by intervening CNM. The association of CNM, CO, and Hαemissions suggests a close link between the colder, denser Hiphase and star formation in NGC 6822. 

Abstract In addition to occupying the extreme, diffuse tail of the dwarf galaxy population, ultra-diffuse galaxies (UDGs) are themselves a key laboratory in which to study star formation in extreme low-density environments. In the second paper of this series, we compare the spatially resolved star formation activity of 22 H i -selected UDGs and 21 “normal” dwarf galaxies within 120 Mpc to predictions within the pressure-regulated, feedback-modulated (PRFM) theory of star formation. To do so, we employ a joint spectral energy distribution fitting method that allows us to estimate star formation rate and stellar mass surface density from UV-optical imaging. We find that the PRFM framework extends successfully to the UDG regime—although the UDGs in our sample show unusually low star formation rate surface densities given their H i content, this low star formation efficiency can be naturally explained by the diffuse structure of the UDGs. In fact, when cast in the PRFM framework, the relationship between midplane pressure and star formation in the UDG sample is in good agreement not only with the “normal” dwarf reference sample, but also with measurements from more massive galaxies. Our results suggest that despite their low star formation efficiencies, the H i -rich UDGs need not be forming stars in an exotic manner. We also find that the UDGs are likely H 2 poor compared even to the overall dwarf population.